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The non-linear dynamics of bending instability and 
vertical structure of a galactic stellar disc embedded 
into a spherical halo are studied with N-body nu- 
merical modelling. Development of the bending in- 
stability in stellar galactic disc is considered as the 
main factor that increases the disc thickness. Corre- 
lation between the disc vertical scale height and the 
halo-to-disc mass ratio is predicted from the simula- 
tions. The method of assessment of the spherical-to- 
disc mass ratio for edge-on spiral galaxies with a small 
bulge is considered. Modelling of eight edge-on galax- 
ies: NGC 891, NGC 4738, NGC 5170, UGC 6080, 
UGC 7321, UGC 8286, UGC 9422 and UGC 9556 is 
performed. Parameters of stellar discs, dark haloes and 
bulges are estimated. The lower limit of the dark-to- 
luminous mass ratio in our galaxies is of the order of 
1 within the limits of their stellar discs. The dark haloes 
dominate by mass in the galaxies with very thin stellar 
discs (NGC 5170, UGC 7321 and UGC 8286). 

Key words: galaxies: haloes - galaxies: spiral - galax- 
ies: structure - dark matter 



o 
u 

6 

u 



> 

m 

\o 

O 

o 



X 



Astron. Nachr./AN 32X (200X) X, XXX-XXX 



Numerical modelling of the vertical structure and dark halo 
parameters in disc galaxies 

A. Khoperskov\ D. Bizyaev^'^, N. Tiurlna^, M. Butenko' 

(1) Volgograd State University, Volgograd, 400062, Russia 

(2) Sternberg Astronomical Institute, Moscow, 119992, Russia 

(3) Apache Point Observatory and NMSU, Sunspot, NM,USA 



Accepted ???? ???? ??. Received ???? ???? ??; in original form 



7777 7777 77 



Abstract. 



1. Introduction 

The luminous matter reveals only a part of the total 
mass in galaxies. Additional fraction of non-luminous 
matter is required to explain rotation curves (RC 
hereafter) of spiral galaxies. Fitting of galactic ro- 
tation curves does not strictly require a lot of dark 
matter Obviously, a modest dark halo with 1/4 of 
the disc's mass is enough to explain a flat RC at the 
radii range r ^ (2 - 4) ■ L dvan Albada et al. 19"85T l. 
where L is the disc scale length. Numerical mod- 
elling of disc galaxies without a massive spheroidal 
subsystem shows that their stellar discs heat up sig- 
nificantly during the evolution, and finally achieve the 
equiUbrium state with too high velocity dispersion 
dOstriker & Peebles 1973; 'Carlber g & Sellwood 1985] 
lAthanasoula & Sellwood 1986 , 
IBottema & Gerritsen 19971 IFuchs & von Linden 19981 



Khoperskov et al. 2003 1. To explain the low ratio 



of the stellar velocity dispersion to the circular ve- 
locity c^'^yVc, which is observed in many galaxies 
at r > 2L, a spherical subsystem as massive as the 
the stellar disc or even more has to be incorpo- 
rated into the models. More additional arguments 
for the massive spherical haloes come from obser- 
vations of galactic polar rings and from long flat 
RCs observed far outside the limits of stellar discs 
dWhitmo re et al. 1987', 'Reshetn ikov & Sotnikova 20001 
llodice et al. 2003 Combes 2002) 1. 

Mass of the spherical subsystem in a galaxy can 
be determined from its rotation with quite high level 
of uncertainty. Incorporation of radial distribution of 
the stellar velocity dispersion into the modelling helps 



to decrease the ambiguities and allows to constrain 
better the range of parameters of galactic subsystems 
(IBottema 19931 IBottema & Gerritsen 19971 IFux 19971 



Khoperskov et al. 2001 ). Unfortunately, observations of 



the steflar velocity dispersion in galactic discs require 
a sensitive spectroscopy and long integration time, so 
the dispersion still can be estimated for a limited num- 
ber of galaxies. An alternative method of including the 
velocity dispersion into the numerical modelling is to 
utilize the disc thickness that, in a contrary, is easy to 
observe in many edge-on galaxies. Relation between 
the vertical scale height, local disc surface density and 
spheroidal subsystem mass M, enables us to include the 
disc thickness into the modelling. At the same time, if 
the velocity dispersion is close to that required for the 
marginal disc stability, the thickness of the stellar disc is 
tightly connected with the spherical-to-disc components 
mass ratio (IZasovetal. 19911 IMikhailova et al. 20011 
I Sotnikova & Rodionov 2006l l. 



One of important problems in the physics of 
galaxies is the identification of the main factor that 
increases the random velocity from 6 -^ 10 km/s 
(young population) to the values of ~ 20 -^ 150 km/s 
(observed in old stellar population). Besides such 
factors of heating as the scattering of stars by giant 
molecular clouds ( Spitzer & Schwarzschild 1951| 
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Spitzer & Schwarzschild 1953 1 and in the spiral 
arms (Jenkins & Binney 19901, the bending in- 
stability is the most important mechanism of the 
vertical velocity dispersion c^ growth in disc galaxies 
dPoliachen ko & Shukhman 19771 IZasov et al. 199T\ 
(Mikhailovaetal. 20011 ISotnikova & Rodionov 20"03l 
ISotnikova & Rodionov 2"005T l. 
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The bending instability in galactic discs 
has been considered with the help of N-body 
simulations. Essential results were obtained 
when the tidal interaction was taken into ac- 
([Hernquist et al. 1993 



count 



Velazquez & White 1999 
BaiUn & Steinmetz 2003t . 



[Weinberg 1998 



Mayer etal. 2001 



Tidal forces are 



erated by the disc bending and increase the disc 
thickness. This is in agreement with observations of the 
probability of bending and warps in isolated galaxies 
and group objects dReshetnikov & Combes 2002l i. 
The emergence of the bending instabilities was 
investigated by IRahaetal. 19911 [Sellwood 1996', 
ISellwood & Merritt 19941 iPatsis et al. 2002., 



IGriv et al. 20021 [Binney et al.T998l A model of 



the galaxy NGC 2146 that involved the bending insta- 
bility was considered by Griv et al. 1998i The bending 
instability reveals itself on much longer time scale, 
typically 10^ years (ISotnikova & Rodionov 2003 ^ in 
the comparison with the gravitational instability in the 
disc plane. iHunter & Toomre 1969l explained the global 
warps observed in discs of some isolated galaxies by 
bending instabilities. 

The disc scale height in galaxies depends on the 
vertical velocity dispersion c,(r). The latter is depen- 
dent of the radial velocity dispersion Cr{r). A dis- 
cussion about the relation between Cr and Cj was 
opened firstly bv IToomre~1 966 Considering a simpli- 
fied model of infinitely thin uniform self-gravitating 
layer, dPoUachenko & Shukhman 19771 ) found an essen- 
tial condition for a stellar system stability against the 
small-scale bending perturbations: c^/cr > 0.37. The 
dynamics of the bending instabilities was also consid- 
ered by lAraki 19861 where a non-uniform volume den- 
sity distribution in the z-direction was taken into ac- 
count. The conclusion about a lower value of the crit- 
ical ratio c^/cr ^ 0.3 was made in that paper (see also 
discussion in lMerrit & Sellwood 19941 ). 

In this paper we point our special attention to the ra- 
dial distribution of the disc thickness and ratio c^/cr that 
are necessary to provide the disc stability against differ- 
ent kinds of bending perturbations. We mostly consider 
isolated late type galaxies without a prominent bulge. 
The presence of bulge would complicate the analysis 
of bending instability. If there is no bulge in a galaxy, 
one can restore the internal part of rotation curve as- 
suming that it is defined by the disc component only in 



the inner (r :$ 2L) regions ( [Zasov & Khoperskov 2003 1. 
However, we also develop our modelling for two well- 
studied disc galaxies with significant bulge component 
because we can incorporate observations of of the stel- 
lar velocity dispersion for these galaxies. 

In § 2 we describe our N-body model. Evolution and 
features of the bending perturbations are considered in 
§ 3. In § 4 we model structural parameters and rotation 
curves of eight edge-on galaxies and evaluate their best 
fitting parameters focusing at dark halo-to-disc mass ra- 
tio. 



2. Dynamical modelling of stellar discs in 
galaxies 

Our N-body modelling is based on the numerical in- 
tegration of the motion equations for A^ gravitation- 
ally interacting particles as it was described in detail 



gen- by Khoperskov et al. 2003' This system of collisionless 



particles has a shape of a disc (and, optionally, bulge) 
embedded into a spherical dark halo. The halo and bulge 
are represented by fixed potentials. The steady state dis- 
tribution of mass in the bulge £)**' and halo p**' are de- 
fined as: 



r^(h,b) 



(r) 






(1) 



where ^ = Vr^+z^^ and (r, z) are the radial and vertical 
coordinates, k equals to 3/2 for the bulge and to 1 for the 
halo. The dimensional spatial scales for the bulge and 
halo are denoted as ai, and a/,, respectively. Our model 
bulge is defined within a sphere with radius ^ < r™''. 

The initial vertical equilibrium in the disc is defined 
by the Poisson equation 

d0rdr (rd(i>0dr) + d^Q>0dz^ = AnG [g + g'' + g'') (2) 

and by a balance of forces in the vertical direction (in 
the first-order approximation) 

cldg0dz = -gd<l)0dz. (3) 

Here O is the gravitational potential, g is the disc vol- 
ume density and c, is the vertical component of ve- 
locity dispersion. The radial component of Jeans equa- 
tion defines the rotational velocity in the stellar disc 
dValluri 19941 ): 



V^ = 



V^+c^l- c\0c]. + r0gcjd(gcj)0dr + r0c^^d{uw)0dz] , (4) 



where the brackets (. . .) denote the mean value. The last 
term in the figured brackets in (|4]) is the chaotic part of 
the radial u and vertical w velocity components. 

The system of equations (|2])-(|3]) can be reduced to 
the equation for dimensionless disc density f(z; r) - 
q{z\ r)/g(z = 0; r), see lBahcall 19841 

cFf0dz^ + 2dlnc.0dzdf0dz - 10/ (df0dzf + 4nG0c^.o-02zo x 



</(/ + 2efVzo) = 0. 



(5) 



.<■// 



where the surface density is cr = 2g(z - 0) ■ zo, P/, 
gi, - l04-KGrdV^0dr and the disc scale height is zo - 
L fdz- For the case of c- = const and g^f^ - the 
volume density profile of the disc in its vertical direction 
is: 



Qiz) = g{z = 0) ■ sech (z/zo) , 



(6) 



where the vertical disc scale for this case is zo = 
cl jnGcr. 

At the flat part of RC, the last term of equation 
(O decreases the vertical scale zo because of additional 
gravitational potential. Nevertheless, this factor is not 
dominant and the most significant reason for zo decreas- 
ing is that the potential of the halo changes the value of 



C7lcr required for the marginal stability of bending per- 
turbations. 



3. Bending instabilities in stellar discs 



We assume that the disc volume density can be writ- 



ten 



£.(r,z)=e(0,0)/,(r)/,(z), 



(7) 



where the function fr{z) describes the vertical den- 
sity profile. We consider two functional forms for the 
vertical density distribution f^{z) in our modelling: 
exp(-z//z,) and sech^(z/zo), where the scale heights h, 
and zq are functions of time and location in the disc. The 
radial density profile /^(r) is assumed to be exponential: 
fr{r) = exp(-r/L). 

The mass M;, of dark halo was calcu- 
lated inside the maximum disc radius ^max- 
We assume that R^ax ~ 4L in our modelling 
(in accordance to Ivan der K ruit & Searle 1981 , 
Ivan der Kruit & Searle 19821 Pohlen et al. 2002, 



HoUey-Bockelmann & Mihos 200 l|i. We also designate 



the spherical-to-disc mass ratio in the galaxies as ju. 

As the initial templates in our models, we consider 
axisymmetric systems in the state of equilibrium. Initial 
radial distributions of the velocity dispersions in radial 
(cr) and azimuthal (c^ = c,cb/2Q) directions keep the 
templates in gravitationally stable state. We assume that 
the stellar disc population is shaped as a quasi-stationary 
balanced system and that the radial and vertical random 
motions of particles keep the disc in the marginal stabil- 
ity against bending waves and perturbations in the disc 
plane. Then we observe the evolution of the discs during 
10^0 turns at its outer edge. 

We use TREE-code approximation in the modelling 
(described in detail by Khoperskov et al. 2003[ ) with 
A^ = (0.5 - 8)10''. Also, in § 3.1.1 we verify some 
model conclusions with independent direct particle- 
particle approach (with less number of N) as well de- 



scribed by,Khoperskov et al. 2003 



For the comparison between the modelling and ob- 
servations, we incorporate both structural parameters of 
stellar discs and their rotation curves at the same time. 
We take into account that the observed structural param- 
eters in edge-on galaxies are results of integration along 
the line of sight. To obtain the edge-on surface density, 
we integrate the volume density of the model disc along 
the line of sight: 



(t{x, z) = r g{^, z) dy (^ = ^x^+y^) . 



(8) 



-yo 



Here x is the projected distance to the centre of the disc. 

The similar approach is applied for the integration 
of RC along the line of sight (it was described by 



3.1. The global bending instability 

The linear analysis of stability against 
global bending perturbations was discussed 

and 



Zasov & Khoperskov 2003 1. 



the 
by 
Polyachen ko & Shukhman 1979| and generalized by 
Vandervoor t 19911 It was shown that the large scale 
instability modes with the azimuthal numbers m - 1,2 
and 3 generate the instability in coUisionless ellipsoidal 
systems in dependence of their axes ratio, velocity 
dispersion and mass of their spherical subsystems. 
The most important result from the papers mentioned 
above is that the area of stability in the parametric 
space depends on the spherical subsystem mass. Once 
the halo is massive enough, the systems with low 
velocity dispersion can be unstable. This gives an 
example of destabilizing role of the halo. On the other 
hand, the higher mass of the spherical subsystem 
requires the lower ratio of Cr/V™" that is necessary 
to keep the disc gravitationally stable. As a result, the 
destabilizing role of the massive halo related to the 
large scale bending perturbations may not reveal itself 
JZasov et al. 199IIIMikhailova et al. 200TT ). 

Our model stellar discs at the initial moment (t - 0) 
are axisymmetric, in the equilibrium state along both the 
radial and vertical directions and gravitationally stable 
in the disc plane. The latter is made by assigning a high 
value to the radial velocity dispersion c, ^ (1 .5 -^ 4) ■ cj, 
where the Toomre stability parameter Qj = Crjcj is 
a function of r ( [Khoperskov et al. 2003 [ ?). In general, 
our initial disc templates can be either stable or unstable 
against the bending perturbations in dependence of Cj.(r) 
distribution. 

The evolution of systems with a small initial value 
of c^Jcr (i.e. systems that are dynamically cold in the 
vertical direction) reveals a gradual growth of the global 
bending instability, which heats up the disc in the verti- 
cal direction and increases its thickness. 

3.1.1. The axisymmetric bending mode (m = 0) 

Development of unstable axisymmetric bending mode 
m - leads to significant disc heating in the verti- 
cal direction (the same conclusion has been made by 
ISellwood 19961 1. The evolution of the vertical coordi- 
nate of the local centre of mass ^(r, t), vertical veloc- 
ity dispersion c^{r, t), disc scale height zo{r, t), disper- 
sion ratio a~ = c^ir, t)/cr{r, t), radial velocity disper- 
sion Crir, t) and rotational velocity V{r, t) in the stellar 
disc where the axisymmetric mode is developing can be 
seen in Fig. [T] We fix the units hereafter assuming that 
Md = 1, G = 1 and AL - 1. In this case the circular 
velocity Vc ~ 1 .5 at the periphery of the disc and the 
corresponding dimensionless rotation period is r a; 4. 

As it is seen in Fig. [T] parameters of the disc stay 
almost unchanged for the first 2.5 turns (t ^ 10 in our 
units). This time the instability modes are being devel- 
oped at the linear stage of the bending instability evo- 
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lution. After f ^ 10, the linear stage is superseded by 
the non-Unear development of the bending instability, 
and the centre of mass ^(r) oscillates with larger ampli- 
tude in the z-direction (see Fig.[T^). The amplitude ^(r) 
rapidly increases in the central regions of the disc up 
to its maximum value and then decreases down to the 
initial state (see curves 1-5 in Fig.[T^), whereas the in- 
creasing of ^(r) occurs much slower at the outer parts 
of the disc (curves 1 1-19). The rapid growth of the ve- 
locity dispersion c, and scale height zq begins after a 
time delay from ^(r), see Fig.[T}5 and c. A flare of stellar 
disc emerges in its central region and then propagates 
outside, towards the periphery. The value of c^/cr rises 
mostly because of c. growth, and in less power due to 
c,- decreasing (Fig.[TJl and e). The azimuthal component 
of the velocity dispersion c^ follows c^, so the relation 
c^p ^ c,se/2Q. is valid almost everywhere. This equation 
can be failed either in a very thick disc (zo/L ^ 0.4), 
or in the very central region, where c,-/c^ < 2Q.lse, that 
means the anisotropy is lower in those regions because 
of the system spherization. Decreasing of c,- and c^ is 
the result of conversion of the kinetic energy of ran- 
dom motions in the disc plane to the kinetic energy of 
random vertical motions. Note that development of this 
bending mode takes place in the disc which is initially 
axisymmetric over all parameters. 

The value of ^ has opposite signs at the disc centre 
and periphery in the considering model (see Fig. [T^). 
This generates a "mexican hat-like" structure as a result 
of this instability (Fig.|2|. The isolines of ^ have concen- 
tric shape. The intermediate regions and periphery at the 
moment t - 25 A ^ 6t are shifted off the disc plane to- 
ward the opposite directions, whereas the central region 
(r $ 0.4) has returned to its initial position, see Fig.|2}5. 
The disc vertical structure driven by the bending mode 
m = is shown in Fig. |2]l. We can see a formation of 
box-shaped features in the disc, which is typical for all 
moments of evolution when the axisymmetric bending 
mode dominates. Note that this feature is not relevant to 
a bar. The considered instability may be responsible for 
the formation of central thickening of discs (i.e. bulges), 
and for the structure of SO galaxies (see Fig.|2ji)). 

A gradual change of the velocity dispersions c, and 
Cr and rotational velocity V begins when the amplitude 
of the vertical oscillations ^ grows significantly. The 
vertical disc heating is accompanied by the increasing 
of c^jcr, but favourable conditions for developing the 
bending instability disappear at certain values of c-jcr- 
As a result, new stable and thick disc forms. The specific 
time for this process depends significantly on the model 
parameters and is of the order of ten turns of the disc's 
outer regions. If the initial distribution of a, = c^jcy 
was subcritical, the linear stage takes longer time due 
to low increment of instability, and because of the heat- 
ing at the non-linear stage was feeble. The amplitude of 
the bending mode, and even the ability of its emergence, 
depends significantly on the initial radial distribution of 



An important consequence is that the final distribu- 
tion of a^{r), as a result of bending instability, depends 
on its initial value a~{r, t - 0). After development of the 
global bending instability, those thin discs that had a low 
value of the scale height ratio gain larger ratio c^ /Cr than 
it would required to keep the discs stable. The explana- 
tion can be found in the essentially non-linear nature of 
the disc heating. As a result, the parameter c./c,- "jumps 
over" the "border of stability" and comes to relaxation 
well above its marginal threshold. 

We verified the numerical effects of our TREE- 
code by reproducing the most important numerical ex- 
periments with a direct particle-particle (PP hereafter) 
code, see description by [Khoperskov et al. 2003| The 
analogue of Fig. [T]that was obtained with the help of PP 
code that started from the same initial conditions and 
was evaluated with only 5 x 10"* particles is shown in 
Fig. [3] The residuals in Fig. [3] are greater than those in 
Fig.[T]due to less number of particles. Nevertheless, all 
features from Fig.[T]can be seen in Fig. [3] 

If the model disc is not massive enough to provide 
the gravitational stability against the bar-mode and mass 
of the halo is low (fj. <l), the main reason for heating 
the initially thin cold disc is the axisymmetric bending 
mode (m = 0) whereas the modes m - 1,2 do not 
emerge at all. Note that the mode m - Q may emerge 
far from the centre (r ; L) and develop there without 
penetrating into the central part of the disc. It can hap- 
pen because of two reasons: either the initial disc had 
Czlcr > ai"' (see § 3.4) and was stable in the central 
regions whereas periphery of the disc is thin and un- 
stable (Fig. m, or the galaxy harbours a concentrated 
and massive bulge. In such models the bending modes 
in the centre have low amplitude. The presence of a 
bulge and/or a halo plays a stabilizing role. Hence, hav- 
ing all other conditions equal, the discs of galaxies pos- 
sessing bulges have to look thinner Note that since we 
did not consider dynamically non-stationary bulges in 
our model, this our conclusion requires further investi- 
gation. 

3.1.2. Bending modes m- 1,2 

In the models with massive halo (/i ^ 2), the axisym- 
metric bending modes m - 1 and m - 2 may develop 
themselves and heat up the stellar disc along its verti- 
cal direction. The surface density of the disc remains 
axisymmetric. The evolution of initially thin disc for a 
case of a heavy halo // = 4 is shown in Fig. |5] In this case 
the bending mode m = 2 is developed initially in the in- 
ner disc regions (r :^ 2L), see the saddle-like features in 
Fig.|6] Fig.|6]shows the distribution of ^ in the disc plane 
for certain moments between t - Q and f = 40 (like in 
Fig. |5]l. The vertical heating and disc's flare due to the 
mode ;« = 2 are very modest in this case. However, a 
non-linear one-arm asymmetric mode m - 1 emerges 
after f ^ 3 and the vertical heating gets more signifi- 
cant. The third stage of the heating starts at f ^ 10 when 



the mode m - begins to dominate all over the disc. 
The growth rate of the vertical dispersion c,- is espe- 
cially high at this time: the disc scale height increases 
by a factor of 2 - 3. The example considered shows the 
process of transformation between the bending modes 
and transition from the asymmetric mode m = 2 to the 
axisymmetric one. In the case of a low-massive spher- 
ical subsystem of low mass, the favourable conditions 
for emerging the global bending instability get worse 
and the disc warp does not emerge in the disc for at least 
20 rotation turns if the initial c- = c-/cr is high enough. 
If the initial template was too far from the stability 
limit, the disc "jumps over" the stability limit, and then 
its final state essentially exceeds this limit. Thus finally 
produces overstabilized and too thick disc. Since such 
stellar discs are not obviously observed, we have to as- 
sume that the initial state of the discs was not too far 
from the stability limit. Note that it agrees with a grad- 
ual building of stellar discs from gas or via minor merg- 
ing when the initial discs do not emerge instantly but are 
built for a continuous time from small fragments, which 
adjust to the discs stability limit. 

3.2. Bending of bars 

In the course of experiments with low massive halo 
(/z < 1 .5), if the initial state of disc is gravitationally un- 
stable, a bar can emerge. The bar formation is accompa- 
nied by its warping (as reported bv' Rahaetal. 19911 1. 

Bending of the bar can be a result of global instabil- 
ity of the bar-mode during its initial stage when the bar 
forms in the initially thin cold disc. Fig.|2]shows the bar 
bendings that increase the vertical velocity dispersion. 
Amplitude of the bar warps decreases essentially when 
the bar thickness increases. We stress that once the bar 
has been formed, it stops the further possible developing 
of the global bending modes, and first of aU it destroys 
the axisymmetric instability with mode m = 0. 

The bar formation runs faster while the initial disc 
is cold (i.e. for small values of Toomre's parameter Qt) 
and as a result, amplitude of bar's bending mode in- 
creases (Fig.|7]l. If the initial disc lays in a slightly sub- 
critical state (cr is just below the level that provides the 
stability to the global bar-mode), the bar emerges very 
slowly, and it is stable against the bending perturbations, 
i.e. the bending of a bar does not necessarily emerge in 
the modelling. 

3.3. The ratio cjcr 

The key parameter responsible for stability of the stel- 
lar discs is the ratio a, = c,/c,-. To stabilize the 
global bending instabilities in the case of low val- 
ues of relative mass of the spherical subsystem in 
a model, the value of c,/c,- has to be greater than 
0.3 - 0.37, as it was figured out from linear analy- 
sis of simple models ( IPoliachenko & Shukhman 19771 
lAraki 1986|[Merrit & Sellwood 1994l i. 



Let's consider a bulgeless model with a moderate 
halo fj. = I, see Fig. [1] and |2] The initial (curve 1) 
and final (curve 2) distributions of a~{r) are shown in 
Fig. [8^. The vertical heating generated by the axisym- 
metric bending mode is as strong that the averaged over 
the disc ratio {c,/cr} = 0.74. 

Once the relative mass of halo is chosen to be higher, 
the ratio c,/c,- needed to marginally stabilize the bend- 
ing perturbations (we refer to it as the critical ratio af"' 
hereafter) gets less. One can see the initial (curve 1) and 
final (curve 2) distributions of a,{r) for the case of very 
massive halo with fi - 4 in Fig.[8j5. The curve 3 shows 
the initial distribution of a^{r), which is necessary to 
provide stability against the global bending modes. At 
the same time, at the disc periphery 0.27 < c^/cr < 0.37. 

A distinctive feature of the considered models at the 
threshold of bending stability is the non-uniformity of 
c~/cr along the radius (see Fig. [8]). For the case of mod- 
erate halo (// ^ 1) with the initial distribution of Cr(r) 
that suppresses the bar instability, the critical value of 
a"'" is a descending monotonous function of r. Its value 
ranges from 0.5 ^ 0.6 at the central regions to 0.3 H- 0.4 
at the periphery (see Fig.[3]l. 



4. Modelling of selected edge-on galaxies 

In order to compare our model predictions with 
observations, we choose eight spiral edge-on 
galaxies. Structural parameters for five of them: 
NGC 4738, UGC 6080, UGC 8286, UGC 9442 and 



UGC 9556, were taken from [Bizyaev & Kajsin 2004 



Superthin edge-on galaxy UGC 7321 was studied by 
[Matthews 1 9991 [Matthews 20001 Radial distributions of 
the stellar velocity dispersion are available in addition 
to rotation curves and structural parameters for two 
large and nearby galaxies NGC 891 and NGC 5170 
jvan der Kruit & Searle 19"8T1 [Morrison et al. 19971 
IBottema et al. 19911 IBottema et al. 1987l i. All the 
galaxies except NGC 891 and NGC 5170 indicate 
a negligible presence of bulge, which simplifies the 
modelling. 

We evaluate the following free parameters of the 
model: yU = MJMd, the radial scale of the halo a/, and 
the disc central surface density o-q. We look for the op- 
timal agreement between the calculated and observed 
disc thickness and rotation curve simultaneously. Here 
Md is the galactic disc mass; Ms = Mi, + Mi, is the mass 
of the spherical component which comprises of the halo 
and bulge in general case. Note that since almost all our 
galaxies have no visible bulge, their spherical compo- 
nent in the most cases means the galactic dark halo. 

We assume that the density distribution in stellar 
discs and bulges (if used) follows the brightness distri- 
bution. It corresponds to assumption of constant mass- 
to-light ratio. 

The structural parameters of the galaxies are figured 
out for the Hubble constant Hq - 75 km s ' Mpc~' . The 
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Table 1. Observational parameters of selected edge-on 5.2 ■ 10^" Mq, cto = 12OOM0/pc^ and Vc,rf/Vcl/=2.2L 
galaxies 0.87. 



Name 



D 

Mpc 



L 
kpc 



(Zo) 

kpc 



kpc 



kpc 



Taf-l-l. UGC 6080 



NGC 891 
NGC 4738 
NGC 5170 
UGC 6080 
UGC 7321 
UGC 8286 
UGC 9422 
UGC 9556 



9.5 

63.6 

20.0 

32.3 

10.0 

4.8 

45.6 

30.6 



4.9 
4.7 
6.8 
2.9 
2.1 
2.0 
3.5 
1.5 (3.6)" 



0.98 
1.30 
0.82 
0.69 
0.17* 
0.26 
0.80 
0.51 



0.49 
0.7 

0.48 
0.14'' 
0.13 
0.51 



21.0 
19.2 
26.2 

9.9 
8.15 

8.0 
14.6 

9.0 



-jThe 
2 UGC 



Vc and V'" in the left panels of Fig. |9] for 
6080 are designated by the solid and dashed 



" two exponential discs are considered: HSB and LSB; the 
latter one has larger scale length. 

* the scale is given for the disc's periphery in case sech(z/z^./,). 
'^ the scale is given for the disc's centre. 
Here D is the distance to the galaxy, L is the exponen- 
tial disc scale length, {zo> is the scale height for the sech^ 
shape of the vertical luminosity profile, (/?-) is the scale 
height for the case of exponential vertical profile, iJ^ax is 
the disc cut-off radius. References to the rotation curves: 1 - 
ISancisi 1979l |Begeman et al. 1991| 2 -JGiovanelli et al. 19971 
3 - Bottema et al. 1987 4 - Karachents ev & Zhou 19"9n 5 - 
[Matthews 1999. 6 - .Bottema et al. 19861 



name of galaxy, adopted distance D, scale length L, ob- 
served mean scale heights (zo) and {h,} (corresponding 
to sech^ and exp luminosity distributions in the vertical 
direction, respectively), stellar disc radius R^ilx, and ref- 
erence to the source of the rotation curve are shown in 
Tabled 

4.1. Notes on individual objects 

Below we present some results and notices on mod- 
elUng of the galaxies. All masses of galactic compo- 
nents given below are estimated within the limits of 



4.1.1. UGC 9422 

Fig. |9] shows the rotation curves y*' together with the 
model circular velocity curves Vc- V*' are shown for 
// = 0.35 and fi - 0.6 (dashed and dash-dotted curves, 
respectively). The "edge-on" values of Vc are plotted for 
the same yU = 0.35 and fi - 0.6 (dotted and solid curves, 
respectively). The model radial distributions of zo are 
shown in the top right panel in Fig.|9] Model with low 
massive halo (fj. - 0.4, the solid curve) produces too 
thick stellar discs. It rules out the maximum disc model 
in which the central surface density ctq - 1440Mo/pc^ 
and the disc contributes 94% into the circular velocity 
{Vc""^/Vc) at r=2.2L. The acceptable agreement with the 
observed disc thickness can be achieved at ju ^ 0.6 (the 
dotted curve). Haloes with yU ^ 1 make the disc to be 
too thin (the model with fi = 1 .6 is shown by the dashed 
curve in the corresponding panel of Fig. |9]l. The best- 
fitting model has parameters: Mj = 8.6 ■ 10'" Mq, M/, - 



3 curves, respectively. The model disc thickness is shown 

4 in the right panel for /v = 0.36 (solid curve), ju = 0.57 

5 (dotted curve) and fi = 0.9 (dashed curve). The stellar 

6 disc is thicker than it is observed for all models with 
1^ < 0.6. The best-fitting model gives the upper limit for 

Ajhe disc mass Md = 5.3- 10'"Mo and cr,, = 1200 M^/pc^ 
for its central surface density. 



4.1.3. UGC 8286 

Left panel for UGC 8286 shows only V*' for two mod- 
els: /x = 0.7 (solid curve) and /i = 1.6 (dotted curve). 
Right panels show radial distribution of observed disc 
thickness for the same cases (solid and dotted curves, re- 
spectively). It is a good illustration that a rotation curve 
can be successfully explained with a wide range of /u 
whereas the disc thickness (right panel in Fig.|9l) helps 
to choose the proper model. The best-fitting disc central 
surface density is the lowest in our sample that reveals 
the LSB nature of this galaxy. 

4.1.4. UGC 9556 

The radial photometric profile of UGC 9556 cannot be 
described by a single exponent. We find that It is better 
fit by two exponents with the scale lengths 1.5 and 3.6 
kpc ( [Bizyaev & Kajsin 2004[ ). The notation in the left 
panel is kept the same as for UGC 6080. In the right 
panel, we show the model thickness for fi - 0.9 (dotted 
curve) andju = 1.9 (solid curve). The best-fitting model 
with yU a: 0.9 produces a bar which can be responsible 
for the deviations from the exponential law in the inner 
regions of the model disc. 

4.1.5. UGC 7321 

Thin stellar disc of UGC 7321 and a lack of bulge indi- 
cates that the galaxy owns a massive dark halo. The no- 
tation in the left panel is kept the same as for UGC 6080. 
In the right panel we show the model thickness for 
yU = 3.4 (dotted curve) and ju = 1.8 (solid curve). The 
best-fitting model withju - 3.4 has Mj = 0.53 • 10'° Mq 
and Ml, = 1.8 ■ IO^'Mq. The stellar disc contiibutes 
^ 60 % into the circular velocity at r = 2.2L. The central 
surface density in this case is only ctq = 22OM0/pc^. 
This highlights the low surface brightness nature of the 
disc inferred by IMatthews 19991 We investigate a case 
of the bar formation in this galaxy, according to the hy- 
pothesis bv lPohlen et al. 2003 1 and make the conclusion 
that it cannot exist in UGC 7321, according to our mod- 
elling. 



4.1.6. NGC 4738 

This galaxy has a bulge according to the photometric 
data, but its luminosity is less than 3% of the disc's that 
enables us to consider bulgeless models for NGC 4738. 
In the right panel we show the model thickness for 
/z = 0.6 (dotted curve) and jj - 0.3 (solid curve). The 
halo with intermediate mass // S 0.5 is suitable for this 
galaxy. The models with less massive halo (fi = 0.3) 
produce too thick discs for this galaxy. 

4.1.7. NGC 5170 

An extended bulge is required to match observed RC 
for this galaxy. On the other hand, the photometry indi- 
cates a presence of very small bulge in this galaxy. The 
observed shape of RC with rather low resolution so a 
disagreement with the model RCs is not unexpected at 
r ^ 2 kpc. Fig.[TO]shows the RC and radial distributions 
of velocity dispersions calculated for the maximum disc 
model and for the model with yu = 1.6. The latter case 
provides a better agreement with observational data for 
both RC and velocity dispersion. Moreover, the maxi- 
mum disc model requires the main disc thickness to be 
1.1 kpc whereas the observations suggest 0.82 kpc. The 
model with /z = 1 .6 produces the projected mean scale 
height 0.8 kpc. In general, the agreement between RC 
and both velocity dispersion (Fig. [TOb ) and stellar disc 
thickness is matched when fi ^ 1 .5. 

4.1.8. NGC 891 



This galaxy has a budge (IBottema et al. 199T1 
IMorrison et al. 1997] l. and we took it into account evalu- 
ating our numerical models. Published values of the ver- 
tical exponential scale height for old stellar disc varies 
from 400 pc to 650 pc dvan der Kruit & Searle 19811 
Bahcall & Kylalis 1985| IShaw & Uilmore 1989! 

Morrison et al. 19971 IXilouris et al. 1999] ). We assume 



the set of structural parameters inferred for NGC 891 
by Ivan der Kruit & Searle 19811 

The observed and model radial distributions of the 
stellar velocity dispersion (Bottema et al. 1991) . stellar 
disc thickness and rotation curves are shown in Fig.fTTI 
The stellar velocity dispersion in the centre of the galaxy 
reaches 160 km/s (IBottema et al. 19911 1 for stars in the 
bulge, and the radial component of the disc velocity dis- 
persion is Cf < 120 km/s. Existence of a molecular ring 
in the galaxy at r ^ 3.5 kpc enables us to assume that 
the galaxy harbours a bar (ISofue & Nakai 1992] l. 

Once the radial distribution of the stellar velocity 
dispersion is taken into account together with condi- 
tion of marginally stable stellar disc, successful model 
should have yu > 1.5. Note that the use of an alterna- 



tive RC from Khoperskov et al. 2001 [ gives /z = 1 .8 with 
other similar parameters. 

Fig.fTTb shows the radial distribution of the disc ex- 
ponential vertical scale height for the set of MhJMj - 
0.4, 0.7 and 1.7. The stellar disc thickness indicates that 



the dark and luminous masses are roughly equal to each 
other within /^max- The value of // inferred from the disc 
thickness is slightly less than that estimated from the ve- 
locity dispersion probably due to additional factors that 
increase the disc thickness in addition to the bending 
instability. 

The disc mass of 1.5 • lO'^M© estimated with 
the shape of rotation curve ( IBottema et al. 199l1 l corre- 
sponds to a very high value of mass- to-light ratio 13.5 
± 4. Our method yields 5 H- 10 for this ratio. This in- 
troduces an argument against the the maximum disc in 
NGC 891. 

Finally, we consider a case of two stellar discs of 
different thickness in the galaxy. This assumption was 
made by IMorrison et al. 19 97 We found that such an 
assumption almost does not change the final best-fitting 
value of /i. 

4.2. Relative thickness of stellar disc versus relative 
mass of spherical subsystem 

A tight connection between the stellar disc thickness 
Zo/L and relative mass of spherical subsystem follows 
from our modelling, see Fig. [12] Each circle in the 
Fig. [12] corresponds to a model galaxy with different 
parameters of the disc and spherical subsystem. The 
figure includes our modelling of the objects studied 
bv lMikhailova et al. 20011 IZasov et al. 1991l in addition 
to results of the present paper. The upper curve corre- 
sponds to mostly bulgeless models. The lower one des- 
ignates the loci of models with noticeable and compact 
bulges. As it can be seen in Fig. [12] the presence of a 
bulge introduces a little scatter to the < zo > /L^MJMd 
dependence that can be easily parameterized in depen- 
dence of the bulge mass. This relation can be used to 
estimate MJMd if Zo/L is known from observations. 

4.3. Results and discussion 

The parameters of best-fitting models for our galaxies 
are summarized in Table 2. In addition to the evaluated 
model parameters, we estimate the dark-to-luminous 
mass ratio Mdark/Mium (it equals to /i if no bulge was 
assumed in the model. In the case of bulge we derive 
the ratio MjarklMium as Mh/(Md + Mb), where Mb de- 
notes mass of the bulge. As one can see from the table, 
the fractions of dark and luminous matter in more than a 
half objects from our sample are roughly the same. The 
exceptions are the thin galaxies NGC 5170, UGC 7321 
and UGC 8286 which have massive dark haloes and low 
surface brightness discs (low ctq in our modelling). The 
same may be referred to the external low surface bright- 
ness disc of UGC 9556. 

In Fig. [13] we compare the relative disc thick- 
ness (zo/L) and its central surface density (here we 
consider only thin external LSB disc in UGC 9556). 
The figure shows a good correlation between the disc 
relative thickness and its central surface density. Our 
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sample is small but much more data from photometry 
can be incorporated. Thus, the disc central surface 
brightness was compared with their thickness by 



Bizyaev & Mitronova 2002 [Bizyaev & Kajsin 2004 
Bizyaev & Mitronova 2009 A good observational cor- 
relation was found: the thicker the disc, the greater its 
central surface brightness. The latter parameter reveals 
the central surface density of disc and corresponds to ctq 
in Fig. [13] It means that the correlation in Fig.fTSJcomes 
from two independent methods: structural studies and 
numerical modelling. 

The relative thickness of the disc is a crucial 
parameter in our modelling, and is relevant to 
the relative mass of the spherical component. The 
relative thickness may span quite a wide range 
from 3 to 12 (see Ivan der Rru it & Searle 1982, 



de Grijs 1997 



Bizyaev & Mitronova 2002 



de Grij s & van der Kruit 1996, 
TBizyaev & Kajsin 2004, 



Bizyaev & Mitronova 2009). According to our conclu- 



sions it means that the relative mass of the spherical 
subsystem varies significantly from galaxy to galaxy. 

All these values are dimensionless whereas the cen- 
tral surface density of stellar disc depends on the am- 
plitude of rotation curve which comes from observa- 
tions. Hence, our model could reproduce a rather thick 
low surface brightness disc (say, without a spherical 
subsystem at all) or, in a contrary, a very dense stel- 
lar disc surrounded by massive halo which makes it 
to be very thin. As it can be seen from Fig. [13] and 
jBizyaev & Mitronova 2002] such kind of galaxies were 
not found among both normal and low surface bright- 
ness spirals. It suggests that some factors beyond our 
model assumptions regulate the connection between the 
disc surface density and dark halo masses. Thus can be 
a result of evolution as well as a consequence of dark- 
luminous matter interaction on the stages of galaxies 
formation. 

Additional factors as density waves, tidal perturba- 
tions or GMO can also be developed in stellar discs and 
amplify the vertical motions in disc. Nevertheless, de- 
veloping the bending perturbations is the most powerful 
mechanism heating up the discs and our model results 
should not underestimate /j. very much. As is was shown 
by IZasov et al. 20041 our approach works well for the 
stellar discs of spiral galaxies Sa-Sd. 

Structural parameters of considered galaxies 
were obtained for mostly old stellar population 



( Bizyaev & Kajsin 2004) , therefore we do not expect 
large variations of mass-to-light ratio in the vertical 
direction. 

The galactic dust layer in edge-on galaxies 
creates a problem of reliability of structural pa- 
rameters. Nevertheless, the parameters for our 
galaxies were obtained using off-plane regions 
(see full discussion in Bizyaev & Kajsin 2004 and 



jBizyaev & Mitronova 2002 1. It should significantly 



decrease the influence of dust to results of our 
modelling. 



Realistic simulations should include "live" halo 
rather that that approximated by a fixed potential. For 
the case of the "live" halo our N-body simulations illus- 
trate an extreme case and the ratio M^IMd found in our 
models shows the lower limit for the halo-to-disc mass 
ratio. We found an ability to run our model upgraded 
with a "live" halo for one of our galaxies, UGC 8286. 
Figure[T4]is the same as the third upper and right panel 
in Figure|9](i.e. its part that concerns the disc thickness 
in UGC 8286) updated with more realistic N-body sim- 
ulations. The asterisks (upper curve), diamonds (next 
lower), and triangles designate the models with ^haio 
- 2-10^, 4-10^ and 6-10^, respectively, where N/,fl;„ is 
the number of bodies in the dark halo in our N-body 
simulations. Note that the curves for N;,a/o - 4-10^ and 
6-10^ are almost identical. As we can see, the more real- 
istic case of the three considered haloes suggests a more 
massive spherical component than that for the fixed po- 
tential, but the difference is of the order of 10%. Thus, 
we encourage to consider our halo-to-disc mass ratio 
given in Table [2] as a lower limit for the real value of 
Mdark/Mium that is nevertheless not far away from it. 

5. Conclusions 

1. Development of the bending instabilities in stellar 
galactic discs is studied with the help of N-body nu- 
merical simulations. The axisymmetric bending mode 
(;« - 0) is found to be the strongest factor which may 
heat up the disc in the z-direction. If no bar is devel- 
oped, the bending modes m - \ and m = 2 are not 
significant contributors to the process of the disc thick- 
ening. The most significant growth of the disc thickness 
occurs at the initial non-linear stage of the bending for- 
mation. Once the bending was destroyed, the vertical 
heating gets less effective. The lifetime of the bending 
mode m = is higher for the larger relative mass of the 
spherical subsystem //. We show that initially very thin 
discs increase their thickness much more rapidly than 
those started from a marginally subcritical state. As a 
result of the global bending instability development, the 
final values of the vertical velocity dispersion cAr) and 
vertical scale height zoif) can be larger than those re- 
quired for the fulfillment of the disc stability condition. 

2. The critical value of the ratio a"" - c^jcy is a 
function of the spherical subsystem parameters. At the 
threshold of stability, the value of ^^"'(r) decreases with 
the distance to the centre. The value of a"" can be twice 
as less at the periphery in the comparison with a^" 
at the centre. The radial trend of (c-jCrY"' can be ap- 
proximated as c,lcr oc exp(-r/LQ,) with very long scale 
length La ^{5 -6)- L. 

3. The average relative vertical disc scale height {zq}/L 
decreases when the relative mass of the halo/v increases. 
This can be utilized to estimate the mass of the spherical 
subsystem for edge-on galaxies. 

4. In the frames of our approach, we conduct N-body 
modelling of eight edge-on galaxies and compare re- 



Table 2. Inferred parameters of the spherical and disc components in the considered edge-on galaxies 



Name 


ah 


Mh 


M, 


<- 


at 


Mt 


M 


Mdark/Mlum 




kpc 


10'" Mo 


10'" Mo 


Me 


kpc 


10'" Mo 






NGC 891 


5.4 


7.13 


10 


715 


1.09 


1.86 


0.9(1.7) 


0.7(1.2) 


NGC 4738 


10.4 


5.5 


11 


870 


- 


- 


0.5 


0.5 


NGC 5170 


8.67 


19.2 


12.8 


491 


2.43 


1.15 


1.6 


1.4 


UGC 6080 


3.53 


3.01 


5.3 


1230 


- 


- 


0.6 


0.6 


UGC 7321 


2.16 


1.8 


0.53 


220 


- 


- 


3.4 


3.4 


UGC 8286 


3.1 


0.76 


0.47 


203 


- 


- 


1.6 


1.6 


UGC 9422 


3.8 


5.2 


8.6 


1200 


- 


- 


0.6 


0.6 


UGC 9556 


3.7 


1.29 


0.57 (0.81) 


715(139) 


- 


- 


1.1 


1.1 



Name of galaxy, scale length of dark halo, mass of halo, mass of disc, disc central surface density, scale length of bulge (if a 
bulge was included), mass of bulge, spherical-to-disc mass ratio fi, and dark-to-luminous mass ratio. 

The numbers in parentheses for UGC 9556 correspond to its outer disc. For NGC 891, the estimates of /i and Mj„ri/M;,„, that 
were obtained from the observed radial velocity dispersion are given in parentheses. 



suits with published rotation curves and surface pho- 
tometry data. For the cases of NGC 5170 and NGC 891 
we additionally incorporate published data on the stellar 
radial velocity dispersion. The relative mass of spheri- 
cal subsystem (the pure dark halo in most cases) is es- 
timated for all the galaxies. We conclude that the mass 
of the dark halo within the optical limits in a half of our 
galaxies is of the order of their disc's mass. The excep- 
tions are a thin LSB galaxy UGC 7321 whose the dark 
halo encompasses about 2/3 of the total overall mass, 
and two more galaxies NGC 5170 and UGC 8286 with 
thin discs. 

5. The central surface density of stellar discs in the 
galaxies of our sample correlates well with the relative 
thickness of the discs. This correlation may be utilized 
to estimate the thickness of stellar discs in non edge-on 
galaxies. 
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Fig. 1. Evolution of the vertical structure parameters of a model stellar disc for the case of yU = 1, a = L when the 
axisymmetric bending mode is being developed: a) the vertical coordinate of the disc's local centre of mass (, b) 
vertical component of the velocity dispersion c^, c) disc scale height zq, d) vertical-to-radial velocity dispersions 
ratio c,/cr, e) Cr, f) rotational velocity of particles in the disc. Different curves are drawn for the following set of 
distances to the centre: rj = 4L ■ (0.05 j - 0.025), where j is shown by the corresponding curve. All the parameters 
are averaged along the azimuthal coordinate. 
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Fig. 2. The shape of a model disc that experiences a development of the bending instability at the moment t - 25.4 
for the model shown in Fig.[T] The isolines denote the distribution ^(r, ip) in the plane of the disc (a) and in the x-y 
plane ((x,y) (b). The bold solid line represents ^ = 0. c) Isolines of logarithm of the surface density are shown at 
the same moments as in the panel (a). The axial symmetry in the surface density distribution can be seen in the 
figure, d) Edge-on view at the disc on the stage of the global axisymmetric mode development. The vertical/radial 
aspect ratio in Figure is increased by the factor of 4. 
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Fig. 3. Results of the same model as in Fig. [T]but evaluated with the help of PP code (see text). 
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Fig. 4. Initial conditions of this model with fi - 4- provide the central parts of disc with stability, and the bending 
mode m -0 develops itself in the outer regions (r >L). a) Initial stages of evolution of f at different distances from 
the centre, b) Distribution of ^ = through the disc. The thick line denotes the isoline f = 0. c) The heating of a 
model disc without a halo at all, in presence of a bulge (M/, = 0.25 M^, b - 0.2 L), when the mode m = is being 
developed. The ratio c/c^ does not rise significantly in regions of the bulge. The notation is kept the same as in 
Fig. [I] 
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Fig. 5. The same as in Fig. [T]but for the model with ^ - A. Distributions of ^ for the internal and external regions 
are shown in separated two top panels. 
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Fig. 6. Distribution of f (x, y) in the disc plane at dif- 
ferent moments of evolution: a) f = 2.7 - the mode 
m = 2 is developing in the central region, b) f = 3.1 
-m = 2— >m=l reorganization, c) f = 3.6, d) f = 4.2 

- the well-developed one-arm mode m = 1 is seen, e) 
f= 10.1-m= 1 — >m = reorganization, d) f = 16.4 

- the well-developed one-arm mode m - Q emerges, e) 
t - 51 - all perturbations in the centre r < L - 0.25 
have almost disappeared. The axisymmetric perturba- 
tions can be seen in the intermediate regions whereas 
the mode m - 2 dominates at the periphery (f, g). h) 
The time dependence of c^Jcr at different radii. 
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Fig. 7. Formation of bar. a) Evolution of the ratio c^Jcr at different radii in the model with bar. b) Distribution of ^ 
in the plane x-y reveals structure of the bending mode at f = 15. The thin solid line designates f = 0. The dotted 
lines show isodenses of the bar Warping of the bar is the reason for rapid vertical heating in the central region 
(lines 1 and 5 at the range of f = 10 H- 20), see text. 
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Fig. 8. Radial distributions of a, = c~Jcr for a) yu = 1 (see Fig.[T]) and b) // = 4 (see Fig.|5]l. The dotted curves 
(1) are for the initial distribution, solid lines (2) designate the final radial distributions and triangles (3) denote the 
critical levels of a-,. 
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Fig. 9, Comparison of our modelling with observed values for UGC 9422, UGC 6080, UGC 8286, UGC 9556, 
UGC 7321 and NGC 4738 (from top to bottom). Left panels: radial distributions of V,. (solid curves) and V'"', 
i.e. observed rotation curves in edge-on galaxies (dotted curves). Right panels: radial distributions of the projected 
scale height. Curves correspond to models with different ju. See explanation of individual curves in ^4.11 Filled 
squares designate observed values (RC or disc thickness) in all panels. 
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Fig. 10. Contributions of different galactic subsystems into the Vc for two models of NGC 5170: ji - 0.83 (a) and 
// = 1.6 (b). c) The radial velocity dispersions: 1 - from observations bv lBottemaet al. 19871 2 - from the model 
with // - 0.83 and 3 - with fi - 1.6. 
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Fig. 11. NGC 891: a) Rotation curves from HI observations bv IBottema et al. 19911 (circles) and lSancisi 19791 
(triangles) and the velocity dispersion from iBottemaet al. 19911 The model distributions of Vc (J) and c^ (2) are 
shown for the case of MhJMd - 0.46, whereas the curves 3 and 4 correspond to Vc and c,-, respectively, for the 
case of Ml, /Md =1.5. b) Observed (crosses) and model (Unes) projected vertical scale height for the models with 
MiJMd = 0.4 (the curve 5), 0.7 (6) and 1.7 (7). 
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Fig. 12. The stellar disc thickness (zo/L) and the 
spherical-to-disc mass ratio MJMd- The circles desig- 
nate models on the stability boundary with different pa- 
rameters of spherical subsystem and initial conditions 
in the modelling. The upper curve corresponds to bul- 
geless models whereas the lower one shows the models 
with noticeable bulges. 
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Fig. 13. The stellar disc thickness (zq/L) and the disc 
central surface density (in Mq pc^^) for the galaxies 
from our sample (circles). The cross designates our 



Galaxy ( Khoperskov & Tyurina 2001 1 and the solid line 
is the least-square linear regression. 
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Fig. 14. Results of "live" halo modelling for one of our 
galaxies, UGC 8286. This figure is the same as the third 
upper and right panel in Figure |9] with addition of three 
models which are designated by the solid curves with 
asterisks, diamonds, and triangles, which correspond to 
2-10^, 4-10^ and 6-10^ bodies in the "live" halo, respec- 
tively. 



